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Abstract 

Understanding the aging process of aerosol particles is important for assessing 
their chemical reactivity, cloud condensation nuclei activity, radiative properties 
and health impacts. In this study we investigate the aging of black carbon con- 
taining particles in an idealized urban plume using a new approach, the particle- 
resolved aerosol model PartMC-MOSAIC. We present a method to estimate aging 
time-scales using an aging criterion based on cloud condensation nuclei activation. 
The results show a separation into a daytime regime where condensation dominates 
and a nighttime regime where coagulation dominates. For the chosen urban plume 
scenario, depending on the supersaturation threshold, the values for the aging time- 
scales vary between 0.06 hours and 10 hours during the day, and between 6 hours 
and 20 hours during the night. 
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1 Introduction 



Black carbon containing particles, or "soot" particles, are ubiquitous in the 
atmosphere and their role for regional and global climate has been widely rec- 



ognized (IPCC 2007). Since black carbon absorbs light (Horvath and Trier 
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1993), it contributes to the aerosol radiative forcing, potentially partially off- 



igetting the cooling effect of scattering aerosol particles such as sulfates (Menon 



et al. 2002). 



Black carbon containing particles originate from the incomplete combustion of 
carbon containing material, hence emissions from traffic are an important con- 



tributor (Bond et al. , 2004). Other important sources for black carbon include 
biomass burning and the combustion of coal by industrial processes. In this 
paper we focus on black carbon from traffic emissions. Measurements of vehi- 
cle emissions from gasoline and diesel cars show that the emitted particles are 
a complex mixture of many chemical species with the main constituents being 



black carbon and organic carbon (Medalia and Rivin, 1982 Toner et al. , 2006). 



Trace concentrations of ionic and metallic species are also present (Kleeman 



et al. , 2000). The exact composition depends on several factors, including the 



18 fuel type, the operating conditions and the condition of the individual vehicles. 
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During their transport in the atmosphere, the composition of these particle 
emissions are further modified. Coagulation, condensation and photochemistry 



are contributing processes, collectively known as aging (Weingartner et al 



1997). During this aging process the composition of the individual particles 



or, in other words, their mixing states change (Furutani et al. , 2008). This 
impacts the particles' physico-chemical properties including their chemical 
reactivity, radiative properties and health impacts. In particular, the aging 
process can change the particles' hygroscopicity from initially hydrophobic to 
more hydrophylic, and hence change their ability to become cloud condensa- 



et al. 2008 



Cubison et al. 2008 ) 



ijbion n uclei (McFiggans et al. 2006 McMurry and Stolzenburg 1989 Moffet 



This is important as models and observations suggest that wet deposition 
represents 70-85% of the tropospheric sink for carbonaceous aerosol mass 



(Poschl, 2005). As a consequence, to assess the budget and impact of black 



carbon, models need to capture the aging process adequately. Many global 
models have simulated both (fresh) hydrophobic black carbon and (aged) hy- 
drophilic black carbon, which can be considered a minimal representation of 



the black carbon mixing state (Cooke et al. 1999 Lohmann et al. 1999; Koch 



2001 Croft et al. 2005). In such a framework only the hydrophilic black car- 



bon is subject to in-cloud scavenging. The conversion from hydrophobic to 
hydrophilic is frequently modeled as a first-order system with the single pa- 
rameter of aging rate or its inverse, the aging time-scale r which represents the 
time-scale on which a population of black carbon containing particles transfers 
from the "fresh" category to the "aged" category. 

While conceptually simple, the actual value of the aging time-scale r is not 



well constrained. Koch (2001) and Croft et al. (2005) compared different ag- 



ing parameterizations in global models and concluded that the model results 
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critically depended on the respective formulation. Riemer et al. (2004) used 



mesoscale simulations to determine r. They derived the aging time-scale for 
black carbon particles as a function of height and time of the day, which sug- 
gested that assuming a single parameter for the black carbon aging time-scale 
is an oversimplification that will incorrectly estimate the black carbon burden. 
However, even though this treatment allowed more detailed insight into the 
aging process, it was still based on ad hoc aging rules inherent to the modal 



model framework that was used (Riemer et al. , 2003). 



Recently, Riemer et al. (2009) developed a particle-resolved aerosol model. 



PartMC-MOSAIC, which explicitly resolves the composition of individual par- 
ticles in a given population of different types of aerosol particles, so that no 
ad hoc aging criteria needs to be invoked. They applied PartMC-MOSAIC 
in a Lagrangian box-model framework to an idealized urban plume scenario 
to study the evolution of urban aerosols due to coagulation and condensation 
over the course of 24 hours. 



In this study, we build upon Riemer et al. (2009) and present a method for es- 



timating aging time-scales of black carbon containing particles using PartMC- 
MOSAIC, based on the idealized urban plume scenario. We take the particle 



population simulated in Riemer et al. (2009) and use a CCN-based aging cri- 



teria to determine whether each individual particle is fresh or aged at every 
timestep, and how many particles transfer between the fresh and aged cate- 
gories during each timestep. By fitting these results to a first-order bulk model 
of aerosol aging we are able to determine the aging timescale without mak- 
ing any a priori assumptions about the aging process. To our knowledge it 
is the first time that a method is presented for explicitly calculating aging 
time-scales. 

Section [2] introduces the model system. In Section [3] we describe the idealized 
plume scenario that served as a basis for the time-scale estimation. Section |4] 
presents our method for deriving time-scales from our model and Section [5] 
shows the results. We summarize our findings in Section |6l 



76 2 Model description 



77 PartMC-MOSAIC is a particle-resolved model that simulates the evolution of 

78 individual aerosol particles and trace gases in a single parcel (or volume) of 

79 air moving along a specified trajectory. For each particle the mass of each con- 
so stituent species is tracked, but the particle position in space is not simulated, 

81 making this a zero-dimensional or box model. In addition to coagulation and 

82 aerosol- and gas-phase chemistry, the model includes prescribed emissions of 

83 aerosols and gases, and mixing of the parcel with background air. The simu- 
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lation results shown here use around 100,000 particles in a volume of around 
16 cm^ (the precise values vary over the course of the simulation). We regard 
this volume as being representative of a much larger air parcel. The model 
accurately predicts both number and mass size distributions and is therefore 
suited for applications where either quantity is required. Details of PartMC- 
MOSAIC and the urban plume scenario are described in Riemer et al. (2009). 
Here we give a brief summary. 



91 The simulation of the aerosol state proceeds by two mechanisms. First, the 

92 composition of each particle can change as species condense from the gas phase 

93 and evaporate to it. Second, the aerosol population can have particles added 

94 and removed, either by coagulation events between particles, by emissions, 

95 or by dilution. While condensation/evaporation is handled deterministically, 

96 emission, dilution and coagulation are treated with a stochastic approach. 



Coagulation between aerosol particles is simulated in PartMC by generating 
a realization of a Poisson process with a Brownian coagulation kernel. For 
the large number of particles used here it is necessary to employ an efficient 
approximate simulation method. We developed a binned sampling method 
to efficiently sample from the highly multi-scale coagulation kernel (in our 
case the Brownian kernel) in the presence of a very non-uniform particle size 



distribution, which is described in detail in Riemer et al. (2009). 



Particle emissions and dilution with background air are also implemented 
in a stochastic manner. Because we are using a finite number of particles to 
approximate the current aerosol population, we need to add a finite number of 
emitted particles to the volume at each timestep. Over time these finite particle 
samplings should approximate the continuum emission distribution, so the 
samplings at each timestep must be different. As for coagulation, we assume 
that emissions are memoryless, so that emission of each particle is uncorrelated 
with emission of any other particle. Under this assumption the appropriate 
statistics are Poisson distributed, whereby the distribution of finite particles 
is parametrized by the mean emission rate and distribution. 



114 Lastly, we must also obtain a finite sampling of background particles that have 

115 diluted into our computational volume during each timestep. In addition, some 

116 of the particles in our current sample will dilute out of our volume and will 

117 be lost, so this must be sampled as well. Again, we assume that dilution is 

118 memoryless, so that dilution of each particle is uncorrelated with the dilution 

119 of any other particle or itself at other times, and that once a particle dilutes 

120 out it is lost. 



We coupled the stochastic PartMC particle-resolved aerosol model to the de- 
terministic MOSAIC gas- and aerosol-chemistry code (Zaveri et al. , 2008) 
in a time- or operator-splitting fashion (Press et al. 2007 Section 20.3.3). 
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124 MOSAIC treats all the globally important aerosol species including sulfate, 

125 nitrate, chloride, carbonate, ammonium, sodium, calcium, primary organic 

126 aerosol (POA), secondary organic aerosol (SOA), black carbon (BC), and in- 

127 ert inorganic mass. 
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MOSAIC consists of four computationally efficient modules: 1) the gas-phase 



photochemical mechanism CBM-Z (Zaveri and Peters 1999); 2) the Multicom- 



ponent Taylor Expansion Method (MTEM) for estimating activity coefficients 



of electrolytes and ions in aqueous solutions (Zaveri et al. , 2005b); 3) the Mul- 



ticomponent Equilibrium Solver for Aerosols (MESA) for intra-particle solid- 



liquid partitioning (Zaveri et al. , 2005a); and 4) the Adaptive Step Time-split 
Euler Method (ASTEM) for dynamic gas-particle partitioning over size- and 



composition-resolved aerosol (Zaveri et al. , 2008). The version of MOSAIC 



implemented here also includes a treatment for SOA based on the SORGAM 



scheme (Schell et al. , 2001). 



138 3 Idealized urban plume scenario 
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For our urban plume scenario, we tracked the evolution of gas phase species 
and aerosol particles in a Lagrangian air parcel that initially contained back- 
ground air and was advected over and beyond a large urban described 
in Riemer et al. (2009). The simulation started at 06:00 local standard time 



(LST), and during the advection process, primary trace gases and aerosol par- 
ticles from different sources were emitted into the air parcel for the duration 
of 12 hours. After 18:00 LST, all emissions were switched off, and the evo- 
lution of the air parcel was tracked for another 12 hours. The time series of 
temperature, relative humidity and mixing height are shown in Figure [1} 

Initial gas-phase and aerosol particle concentrations as well as gas phase and 



particle emissions were the same as in Riemer et al. (2009). The gas phase 



emissions varied throughout the emission time interval according to a typical 
diurnal cycle found in polluted urban areas. 

The initial particle distribution, which was identical to the background aerosol 
distribution, was bimodal with Aitken and accumulation modes (Jaenicke 



[T993p . We assumed that it consisted of (NH4)2S04 and POA, as shown in Ta- 
ble [1} We considered three different types of carbonaceous aerosol emissions: 
1) meat cooking aerosol, 2) diesel vehicle emissions, and 3) gasoline vehicle 
emissions. The parameters for the distributions of these three emission cate- 



gories were based on Eldering and Cass (1996), Kittelson et al. (2006a), and 



Kittelson et al. (2006b), respectively. For simplicity in this idealized study, the 



particle emissions strength and their size distribution and composition were 
kept constant with time during the time period of emission. 



5 



168 
169 
170 
171 
172 
173 
174 
175 
176 
177 



Furthermore, we assumed that every particle from a given source had the 
same composition, with the species listed in Table [T| since to date the mixing 
state of particle emissions is still not well quantified. In particular, we assume 
that the diesel and gasoline exhaust particles consist exclusively of POA and 



Rivin 1982; Kleeman et al. 2000) 



J3C, w hich is very nearly the case ( Andreae and Gelencser 2006 ; Medalia and 



Figure |2] shows time series of the bulk aerosol mass concentrations as they 
result from this urban plume scenario. We observed a pronounced production 
of ammonium nitrate, reaching nitrate mass concentration of up to 26 fig m~^ 
and ammonium mass concentration of 10 /xgm~^ in the late afternoon. Sulfate 
mass concentrations increased from 4.1/igm~^ to 6.0/igm~^ due to conden- 
sation of photochemically produced sulfuric acid. POA and BC were directly 
emitted (with a temporally constant rate) and accumulated to 11 /igm~^ and 
4.3yUgm~^, respectively, until 18:00 LST when the emissions stopped. After 
18:00 LST the mass concentrations declined due to dilution, especially nitrate 
and BC for which the background mass concentration were zero. 



178 3.1 Characterizing mixing state 



To characterize the mixing state and to discuss the composition of a particle, 
we refer to the BC mass fractions as 

WBCdry = (1) 

A'dry 

179 where /iBc is the mass of BC in a given particle and fidry is the total dry mass. 

180 Based on this quantity, we then define a two-dimensional number concen- 

181 tration that is a function of both particle composition and diameter. The 

182 two-dimensional cumulative number distribution N-Bc,dry{w , D) is the number 

183 of particles per volume that have a diameter less than D and a BC mass 

184 fraction of less than w. The top panels in Figure [3] show the corresponding 

185 two-dimensional distributions, normalized with the respective total number 

186 concentrations, after 1 hour and after 24 hours of simulation time. Since even 

187 at the time of emission no particles were pure BC, particles were not present 

188 at w^^Q i^j-y = 100%. Fresh emissions from diesel vehicles (tygQ ^ry ~ 

70%) and 

189 gasoline vehicles (ifBc,dry = 20%) appear as horizontal lines since particles 

190 in one emission category were all emitted with the same composition. At 

191 WBc,dry = 0% all the particles appear that do not contain any BC (i.e. back- 

192 ground particles and particles from meat cooking emissions that have not un- 

193 dergone coagulation with particles containing BC). After 1 hour (07:00 LST) 

194 a small number of particles between these three classes indicate the occur- 

195 rence of coagulation. Comparing this result to the result for the end of the 
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196 simulation, we note that at the end of the simulation particles with diameter 

197 below D = 0.03 fim were heavily depleted due to coagulation. A continuum of 

198 mixing states formed between the extreme mixing states of WBC,dry = 0% and 

199 WBCdry = 70%. 



200 3.2 Calculating CCN activity 
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Given that we track the composition evolution of each individual particle 
throughout the simulation, we can calculate the critical supersaturation Sc 
that the particle needs in order to activate. We use the concept of a dimen- 



j^ionless hygrosco picity parameter suggested by Ghan et al. (2001) or Fetters 



and Kreidenweis (2007). In Fetters and Kreidenweis (2007) this parameter is 



denoted by k, and we adopt their notation for the remainder of the paper. 
This concept has the advantage that results from laboratory measurements 
can be used to quantify the hygroscopicity of complex compounds for which 
K values cannot be calculated in a straightforward manner. The overall k for 
a particle is the volume-weighted average of the k values of the constituent 
species. This requires the assignment of individual k, values for each aerosol 
component in MOSAIC. 
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Fetters and Kreidenweis (2007) compiled a table (Table 1 in their paper) with 
K values for a variety of inorganic and organic species based on recent labora- 
tory measurements or on thermodynamic model calculations. For (NIl4)2S04 
and NH4NO3 they report k values of 0.61 and 0.67, based on calculations 
by Glegg et al. (1998) and measurements by Svenningsson et al. (2006), re- 



spectively. Based on this we assume k = 0.65 for all salts formed from the 
NH|-S04"-N03 system. For all MOSAIG mo del species that rep resent SOA 
we assume k = 0.1, based on measurements by Frenni et al. (2007). Following 
Fetters et al.l ([2006| we assume k = 0.001 for FOA and k = for BG. The 



critical supersaturation Sc for a particle of diameter D and volume-weighted 
hygroscopicity parameter k is then given by 



Sr 



c 



where 



C = 



4^3 



and 



A 



(2) 



(3) 



27 " " R*TpJ 
with being the surface tension of water, the molecular weight of water, 
R* the universal gas constant, the water density, and T the temperature. 

Similarly to the use of WBc.dry above, we can use Sc to define a two-dimensional 
cumulative number distribution Ns{D, Sc) in terms of size and critical super- 
saturation. The bottom panels in Figure[3]show examples of the corresponding 
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218 two-dimensional distributions after 1 hour (left) and after 24 hours (right) of 

219 simulation. While freshly emitted diesel, gasoline and meat cooking particles 

220 differ in their BC and POA mass fractions, they are very similar in their hy- 

221 groscopicity with initial n values close to zero. After 1 hour they are visible 

222 as the dark band of high number concentrations at high Sc values. Separated 

223 from this we see another dark band representing the most hygroscopic parti- 

224 cles, consisting of wet background particles. They contain the largest fraction 

225 of inorganic mass (ammonium, sulfate, and nitrate), hence their critical su- 

226 persaturation is lowest at a given size compared to the other particle classes. 



227 Directly slightly above the most hygroscopic band at 1 hour is a weaker band, 

228 which represents the dry background particles. Because they are dry and the 

229 vapor pressure of HNO3 is still low, nitrate formation does not occur on these 

230 particles. The coexistence of wet and dry particles can be explained by the 

231 fact that the relative humidity falls below 85% at 06:42 LST, which is the del- 

232 iquescence point of the inorganic mixture of ammonium, sulfate, and nitrate. 

233 Particles that exist before that time contain water and take up nitrate. They 

234 stay wet throughout the whole day as a result of the hysteresis of particle 

235 deliquescence and crystallization. Particles that are emitted after 06:42 LST 

236 do not contain water and take up nitrate only much later. Hence after 1 hour 

237 of simulation, at a given size the fraction of highly hygroscopic inorganics is 

238 higher for the wet particles, which results in a higher n value and lower critical 

239 supersaturation Sc- 



240 The individual bands are not completely separated at 1 hour, but the regions 

241 in between have started to fill out. The reason for this is the occurance of 

242 coagulation, which produces particles of intermediate composition and hence 

243 corresponding intermediate Sc values. After 24 hours the population as a whole 

244 has moved to lower critical supersaturations, and the distribution with respect 

245 to Sc has become more continuous. Given a certain size, the critical supersat- 

246 uration ranges over about one order of magnitude. 



247 Figure |4] shows CCN properties as more traditional CCN spectra. This repre- 

248 sentation is the one- dimensional projection of the bottom panels of Figure [3] 

249 onto the critical supersaturation axis, plotted as a cumulative distribution. 

250 The change in CCN properties over the course of 24 hours is obvious. After 

251 1 hour a supersaturation of S* = 1.5% is necessary to activate 50% of the 

252 particles by number. This required supersaturation decreases to S = 0.1% 

253 after 24 hours. In the following section we use the results of this urban plume 

254 scenario as a basis for estimating the aging time-scales. 
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255 4 First-order Models of Aging 



256 In this section we describe tlie first-order model of black carbon aging to which 

257 we fit the particle-resolved data simulated with PartMC-MOSAIC, in order 

258 to determine the aging time-scale. We emphasize that we do not actually 

259 simulate using the first-order models presented in this section. Such first- 

260 order systems are frequently used, however, to model the conversion from 

261 hydrophobic to hydrophilic black carbon with the single parameter of aging 



270 
271 



272 



280 
281 



262 rate or its inverse, the aging time-scale r (Croft et al. , 2005). More specifically, 

263 the aging time-scale represents the time-scale on which particles that would 

264 initially not activate turn into particles that can be activated, given a certain 

265 chosen supersaturation threshold. Budget equations for the fresh and aged 

266 populations can be formulated in terms of either number or mass. A number 

267 based aging time-scale is relevant for aerosol indirect forcing as the cloud 

268 optical properties depend on the cloud droplet number distribution. On the 

269 other hand, a mass aging time-scale is relevant for in-cloud scavenging and wet 
removal of BC mass. In the following we will present results for both number- 
and mass-based aging time-scales. 



At time t, the total number concentration Abc(^) of BC-particles is the sum 

273 of the number concentration of BC-containing fresh particles Nf{t) and the 

274 number concentration of BC-containing aged particles Aa(t). We define anal- 

275 ogously the total BC mass concentration Mbc(^), the BC mass concentration 

276 in fresh BC-containing particles Mf (t), and the BC mass concentration in aged 

277 BC-containing particles Ma(t). The aged and fresh populations are separated 

278 by applying an aging criterion, in our case activation at a certain supersat- 

279 uration threshold Sc- Fresh particles are those with critical supersaturation 
above the threshold value, while aged particles have critical supersaturations 
below the threshold. 



In the PartMC model we explicitly track a finite number of particles in a com- 
putational volume V. The number of fresh and aged BC-containing particles 
in the volume V(tk) at time is denoted by rifitk) and rig^itk), respectively. 
Similarly, m{{tk) and mg^{tk) are respectively the total mass of BC in fresh and 
aged BC-containing particles in V{tk). The number and mass concentrations 
of fresh BC-containing particles in V{tk) are then given by 

N(f\- ""f^^^^ M(f\- ""^^^^^ 



282 The fresh and aged number and mass concentrations can change due to emis- 

283 sion and dilution, while condensation and coagulation can transfer number 
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and mass concentration from the fresh to the aged population and vice versa. 
Changes in number and mass concentrations also occur due to temperature 
and pressure changes. In our model we neglect at present the impact of hetero- 
geneous reactions on the surface of the particles, although studies have shown 
that these also contribute to the aging process (Rudich et al. , 2007). The gain 
and loss terms for the fresh and aged BC-containing populations are given in 
Table d 



291 To express changes in number and mass for coagulation we consider that all 

292 constituent particles are lost during a coagulation event and the product of 

293 coagulation is a gain of a new particle. Table [3] shows the overview of all 

294 possible combinations and the resulting terms for each of those combinations. 

295 For example, let us assume that there are four independent coagulation events 

296 within a single timestep: one event between two fresh BC-containing particles 

297 resulting in an aged particle, two events between fresh and aged BC-containing 

298 particles each resulting in an aged particle, and one event between a fresh 

299 BC-containing particle and a non- BC-containing particle resulting in a fresh 

300 particle. Then we have losses An^^f t^) = 5 and An'^f(tk-i,tk) = 2, 

301 and gains Anf"^^(tk_i,tk) = 1 and An'^°^^(tk-i,tk) = 3. 

302 Note that for each particle pairing in Table [3] there can be two outcomes. For 

303 example, coagulation of a small fresh and large aged particle generally pro- 

304 duces an aged particle, while coagulation of large fresh and small aged gener- 

305 ally produces a fresh particle. When Sc, as calculated in equations ([2])-(|3]), is 

306 used as the criterion for fresh versus aged, it can be shown that coagulation 

307 of two aged particles never produces a fresh particle, and that the coagulation 

308 of two fresh particles with Sc values close to the cutoff value can produce an 

309 aged particle. 

Coagulation can result in a net loss of number but must conserve mass. We 
thus have 



Anr^(tfc_i,tfc) < Anr^l{h-i,h) + A<!!|(4_i,4) (7) 
Amrs(tfc_i, tk) = Am^rf + Am^-^ifc^i, t^) (8) 



310 and similarly for coagulation resulting in fresh particles. 



The following continuous equations describe the evolution of the number and 
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mass concentrations of fresh and aged BC-containing populations: 

dNi{t) 



dt 



^density + N^^it) - iVf'^''""°^(t) 



(9) 



i\Vcond/ 



7Vaging(-() 



dt 



dMi{t) 
dt 



^density + Nr^\t) - N f (t) 

+ N^zi^t) + Nr^it) - - ( iv-t(t) + N:°:f{t) 



(10) 



^density + M^^i^{t) - Aff 



(11) 



^/dc- aging (J-) ^/aging((-) 



(12) 



A;/aging(j) 



Md<:-aging(t) 
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Note that the mass equations ( 11 )- 



12 ) have been simphfied from forms identi- 



cal to the number equations ([9 )-( 10 ) by using equation ^ for the conservation 
of mass during coagulation. Condensation (or rather evaporation) and coag- 
ulation can in principle also produce a transfer from aged to fresh particles. 
This is reflected by the terms denoted as Ar'^<'-'^s^"s(t) and M^^-'^g^°s(t) above. 

The aging time-scales tn and tm for number and mass are then determined 
by the first-order models: 



1 



-iVf(t) 



M'^si"s(t) = -J—Mf{t) 

TM{t) 



(13) 
(14) 



The discrete versions of the balance equations ([9])-(12) are: 

riiitk) - nf(t,_i) = Anr^it,_^,tk) - Anf "*-'^(t,_i, t,) (15) 

+ A<-f (t,_i, tk) + Anr^(tfc_i, tfc) - Anr^f{tk-i,tk) 
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- ^a(tfe-i) = Anr*(t,_i,t,) - Anf (16) 
+ An?™^(t,_i, tk) + Anrs(tfc_i, t^) - A<-s(t,_i, t^) 

- fA<°_!:f(tfc_i,4) + A<!:f(t,_i,t, 



miih) - mf(t,.i) = Amr*(tfc^i,4) - Amf "*^°"(4_i, t^) (17) 

+ f Am™_:f(tfc_i,tfc) + AmlZKtk-i,t, 



- (Am-°,^(t,_i,t,) + Am?!:f(^fc-i,^fc)) 



ma(tfc) - m,(t,_i) = Amr*(tfc-i,4) - Amf "*--(4_i, t,) 

+ fAm?-;f(t,_i,t,) + Am^!:f(^^__^^^^ 



- ( Am:°_:f(t,_i,t,) + Amr_:f(4_i,tfc) 



316 Here equation (|8| for the discrete conservation of mass during coagulation was 

317 used to simplify the mass equations, as in the continuous case. 

By comparing the continuous equations ([9])- ( 12 ) to the discrete equations ( 15 )- 
(18) using the relationships ^ we see that the continuous aging terms can be 
approximated by: 

318 where the use of V at time tk is because the computational volume V is 



319 updated first within each timestep in the PartMC algorithm (|Riemer et al. 

320 



2009 Figure 1). 



From equations ( [l3) ) and (14) and the relationships the aging time-scales 
tn and tm are then approximated by: 

^ ' An-g-g(tfc_i,tfc) ^ ^ 

(tfc -tfc-i)mf(tfc) 
"^^'^•^ ^ Am--(4_„4)- ^''^ 
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For the analysis in Section [s] we additionally define the aging time-scale r§°°^ 
ignoring the impact of coagulation, an average time-scale rN,day during the 
day, and an average time-scale rN,night during the night, given by: 



1 /-ISiOOLST 



lTN,dayJ ^ = 1^ / (Mt)) ^dt (24) 

04:00 LST i 

[Mt)) dt, (25) 



(T"N,msht) ^ 



10 hours J 18:00 LST 

321 and similarly for r^°°<^, rM,day, i"M,night, ^^N^day and so on 



322 5 Results 



In this section we show results for supersaturation thresholds ranging between 
Sc = 0.1% and Sc = 1.0%, as these are typically achieved in updrafts in 



stratus and cumulus clouds (Warner, 1968). To limit the number of figures we 
use Sc = 0.6% as a base case. Figure |5j shows the time series for A^f, A^a, and 
Nbc = Nf + Ni^ for Sc = 0.6%. The total number concentration A^bc of BC- 
containing particles increased until 18:00 LST due to the emission of particles. 
After the emissions stopped, A^bc decreased as a result of continued dilution 
and coagulation. The time series for A^f and A'a show that both increased 
in the morning hours. A pronounced transfer from fresh to aged occurred 
after 11:30 LST, the time when nitrate formation started taking place on the 
dry particles (compare Figure [2]). This process efficiently contributed to the 
conversion of fresh particles to aged particles, which was reflected by a short 
aging time-scale. 

The left column of Figure [6] compares the aging time-scales and r™°'^ 



computed according to equations (21) and (23) for different supersaturation 



thresholds. The top, middle and bottom panels shows the results for super- 
saturation thresholds Sc = 0.1%, Sc = 0.6%, and Sc = 1%, respectively. The 
grey shading is the raw data, and the black lines are results computed from 
smoothing the transfer rates with a Hann filter with a window width of 1 hour. 

The solid lines represent tn, including the contributions due to both coagula- 



tion and condensation according to equation (21). For Sc = 0.6%, tn started 
off in the morning with values of tn ~ 20 hours. It decreased sharply after 
11:00 LST, which is the time when photochemistry was at its peak, and ni- 
trate formation was most pronounced, hence leading to a fast aging process. 
Between 12:00 LST and 16:00 LST, tn was less than 1 hour, and reached val- 
ues as low as 0.2 hours. After 16:00 LST, as photochemistry slowed down, tn 
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349 increased again, reaching a plateau of tn ~ 10 hours during the evening and 

350 night. 



The broken hues represent r^°"*^ as defined in equation (23), i.e. the contribu- 
tion due to coagulation is ignored. There was only a small difference between tn 
and r^°"'^ during midday and early afternoon when condensation was operat- 
ing very effectively. However during morning, afternoon and night, neglecting 
coagulation lead to larger time-scales (up to one order of magnitude). 

A similar pattern was found for Sc = 1%, shown in the bottom panel. For 
Sc = 0.1% we generally obtained larger time-scales. During the morning, tn 
was about 50 hours and decreased to 10 hours in the early afternoon. There 
was a short period around 16:00 LST when tn dropped to 2 hours. Obviously, 
at Sc = 0.1%, even after the growth due to condensation of ammonium nitrate, 
the particles were still too small to be activated to the same extent as seen 
for the larger supersaturations. During the following night tn was around 10- 
20 hours. For this low supersaturation the contrast between day and night was 
not as pronounced as for Sc = 0.6% or = 1%. 

Qualitatively, the temporal evolution of tm (right column) was similar to the 
number-based result tn. However the day /night contrast was more pronounced 
for the case with Sc = 0.1%, and the time-scales based on mass during the 
day were lower than the ones based on number. 

Figure [7] shows the individual transfer terms of BC-containing particles for the 
case Sc = 0.6%. The transfer rate N^'^f due to coagulation from aged to fresh 
was very small throughout the whole day, remaining below 0.01 cm~^s~^. The 
transfer rate A^f^^f due to coagulation from fresh to aged followed the time 
series of Af. The minimum of A^f during the early afternoon was reflected in 
a minimum of the transfer rate N^!^^- The transfer rate N^^^ due to conden- 
sation from fresh to aged was large between 11:30 and 15:00 LST, consistent 
with the decrease of tn in Figured, Lastly, there was a non-zero transfer N^'^f 
from aged to fresh due to condensation, which was larger than N^^^ towards 
the end of the simulation. This was related to a shrinking of the particles due 
to decreasing relative humidity (compare Figure [T]). 

Figure [8] summarizes the results for the different aging time-scale deflnitions. 



381 We calculated rN,day and rN,night according to equation (23 ) for several hundred 

382 different supersaturation thresholds between Sc = 0.1% and Sc = 1%. We 

383 also distinguished between the deflnition of r including the transfer due to 

384 coagulation and condensation (solid lines), and including only the transfer 

385 due to condensation, t^'^^j and T]^°JJ?^ht (broken lines). 

386 For this particular urban plume scenario the following general features emerge: 

387 during the day, condensation of semi- volatile species, in our case especially am- 

388 monium nitrate, was the dominant process for aging. The time-scales based 
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389 on number were larger than the time-scales based on mass during the day, 

390 by roughly a factor of five. The aging time-scales had a strong dependence 

391 on supersaturation threshold. For Sc = 0.1%, rN,day was 10 hours, whereas for 

392 Sc = 1%, TN,day was ouly 0.2 hours. During the night condensation was limited, 

393 hence coagulation was the dominant aging process. For low supersaturation 

394 thresholds, the time-scales based on number were smaller than the time-scales 

395 based on mass, but this difference decreased for higher supersaturation thresh- 

396 olds. 



397 6 Conclusions 



In this paper we presented a method for explicitly calculating aging time- 
scales of black carbon aerosol using particle-resolved model simulations with 
PartMC-MOSAIC. We developed number-based and mass-based aging time- 
scales using the activation of the particles at a given supersaturation as a cri- 



J^erion for aging. We applied this method to an urban plume scenario (Riemer 



et al. , 2009) for a range of supersaturation thresholds between 0.1% and 1% 

404 and considered condensation of secondary substances and coagulation as aging 

405 mechanisms. Aging due to heterogeneous processes (chemical aging) were not 

406 included. 

407 For this particular scenario we found a separation into day and night regimes. 

408 During the day the condensation-induced aging dominated, in particular due 

409 to the formation of ammonium nitrate. Therefore the condensation-only number- 

410 based aging time-scale T^°^^y was almost the same as the total time-scale rN,day 

411 and similarly for mass. The daytime aging number-based time-scale rN,day was 

412 about 10 hours for a supersaturation threshold of Sc = 0.1% and decreased to 

413 0.06 hours for a threshold of Sc = 1%. The daytime mass-based aging time- 

414 scale rM,day was about a factor of five lower than rN,day for all supersaturation 

415 thresholds. 

416 During the night, the absence of condensable species caused the number-based 

417 aging time-scale rN,night to be about one order of magnitude larger than during 

418 the day. Coagulation became dominant, which was reflected by the fact that 

419 the condensation-only time-scale T^°ni^ht "^^^ an order of magnitude larger 

420 than the total time-scale rN,night- The nighttime aging time-scales therefore 

421 depended on the particle number concentrations. We suspect that chemical 

422 aging would have its largest impact during periods when condensation is not 

423 dominant, i.e. during the night in our case. 

424 Compared to the time-scales used in global models, which are typically on the 



order of 30 hours (Chung and Seinfeld, 2002 Koch, 2001, e.g.), our time-scales 



were much shorter, in particular during the day. This confirmed findings by 
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427 Riemer et al. (2004) who showed with a completely different approach that 

428 daytime and nighttime aging regimes exists, and that aging during the day 

429 proceeds very rapidly. 



430 However some caveats need to be emphasized: our urban plume scenario repre- 

431 sents only one scenario, with very polluted conditions and fairly high number 

432 concentrations during the night [N^c ~ 5000 cm~^). For lower number con- 

433 centrations, we expect the aging time-scales during the night to increase. 
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Initial/Background 


N{ 


m 




(1) 


Composition by mass 


Aitken Mode 


3.2 


•109 


0.02 


1.45 


50% (NH4)2S04, 50% POA 


Accumulation Mode 


2.9 


•109 


0.116 


1.65 


50% (NH4)2S04, 50% POA 


Emissions 


E (m 


-2s-l) 


Dgn (Aim) 




Composition by mass 


Meat cooking 


9- 


10^ 


0.086 


1.9 


100% POA 


Diesel vehicles 


1.6 


•108 


0.05 


1.7 


30% POA, 70% BC 


Gasoline vehicles 


5 • 


10^ 


0.05 


1.7 


80% POA, 20% BC 



Table 1 

Initial and emitted aerosol distribution parameters. The initial aerosol distribution 
is also used as the background aerosol distribution. The percentages for the com- 
position are by mass. E is the area source strength of particle emissions. Dividing 
E by the mixing height and multiplying by a normalized composition distribution 
gives the number distribution emission rate. 
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Terms 



Description 



Ni{t), N.,{t) 



^dilution (^^^^ ^dilution^^^ 



^density^^^^ ^density 



Number concentration of fresh/aged BC-containing 
particles. 

Gain rate of number concentration of fresh/aged 
BC-containing particles due to emission. 

Loss rate of number concentration of fresh / aged 
BC-containing particles due to dilution. 

Gain rate of number concentration of aged/fresh 
BC-containing particles due to condensation or 
evaporation on fresh/aged particles. 

Gain rate of number concentration of fresh/aged 
BC-containing particles from coagulation events. 

Loss rate of number concentration of fresh 
BC-containing particles to coagulation events 
resulting in fresh/aged particles. 

Loss rate of number concentration of aged 
BC-containing particles to coagulation events 
resulting in aged/fresh particles. 

Gain rate of number concentration of fresh/aged 
BC-containing particles due to air density changes. 



Net transfer rate of fresh-to- aged/aged-to- fresh 
number concentration of BC-containing particles. 

Table 2 

Description of individual terms in equations ^ and (10). With the exception 



of iV^'^"^^*^(t) and N^'^'^'^^^^ (t) all of these terms must be non-negative. The same 



rdensity / 



notation is for the terms in equations (11) and (12) for mass concentration, and for 



the corresponding discrete equations (15)-(18). The discrete terms are expressed 



as a change in number or mass within a timestep, so that Anf^^(ifc_i, t^) is the 
number of BC-containing particles in the computational volume that change from 
fresh to aged due only to condensation during the timestep from time tk~i to tk, for 
example. There are no discrete terms for air density changes as they are incorporated 
by changing the computational volume V. 
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Particle 1 


Particle 2 


Resulting 
particle 


Non-zero loss terms 


Non-zero 
gain terms 


fresh 


fresh 


fresh 


An^Zf = 2 


An^°^s = 1 


fresh 


fresh 


aged 


An^Zl = 2 


AnT^^ = 1 


aged 


fresh 


fresh 


AnlZI = 1, An™7 = 1 


An^°"^ = 1 


aged 


fresh 


aged 


An^l = 1, An^!:f = 1 


An^"^ = 1 


aged 


aged 


aged 


An^°4| = 2 


Anl'"'^ = 1 


fresh 


non-BC 


fresh 


An^!:f = 1 


An^°^s = 1 


fresh 


non-BC 


aged 


An-^f = 1 


An^^^ = 1 


aged 


non-BC 


fresh 


Ann = 1 


An^°^s = 1 


aged 


non-BC 


aged 


An^-SI = 1 


An^^^ = 1 



Table 3 

The different coagulation events and the resulting expressions for the loss and gain 
terms of the number of fresh and aged BC-containing particles. Similar expressions 
exist for mass changes Am in the particle-resolved model and for the number rates 
N and mass rates M in the continuous model. 
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Fig. 1. Time series of temperature, relative humidity, and mixing height over the 
course of the 24 hour simulation. The pressure and water mixing ratio were kept 
constant. 
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Fig. 2. Time series of mass concentrations of selected aerosol species: nitrate (NO3), 
ammonium (NH4), POA, sulfate (SO4), BC, and SOA. Particle and gas phase emis- 
sions were present from 06:00 to 18:00 LST. 
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Fig. 3. Normalized two-dimensional number distributions after 1 hour (07:00 LST) 
and 24 hours (06:00 LST the next day) of simulation. The top panels show the nor- 
malized value of the two-dimensional distribution d'^N-QQ ,^j.y[D,w)/[d\ogiQD dw) 
with respect to diameter D and BC dry mass fraction t^BCdryi while the 
bottom panels show the normalized value of the two-dimensional distribution 
d'^Na,{D, Sc)/{dlogiQ D dlogiQSc) with respect to diameter D and critical super- 
saturation Sc- 
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Fig. 4. Cloud condensation nuclei spectra after 1 hour (07:00 LST) and 24 hours 
(06:00 LST the next day) of simulation based on the simulation results shown in 
Figure [3} 
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Fig. 5. Time series of total number concentration TVbc of BC-containing particles 
and number concentrations N{ and ATg, for fresh and aged BC-containing particles, 
respectively. The critical supersaturation separating fresh and aged particles is set 
to 5c = 0.6%. 



25 



critical supersaturation St 

■irw, 

\ 



0.1% 



number 

TT — r. 



1 1 





' ''11 




' • ..MM'm 







































critical supersaturation Sc = 0.6% 




critical supersaturation Sc = 1.0% 





critical supersaturation Sc = 0.6% 



J .1 



; y \ TM : 













critical supersaturation Sc = 1.0% 






. Tm ,.; 















06:00 12:00 18:00 00:00 06:00 06:00 12:00 18:00 00:00 06:00 

local standard time (LST) (hours:minutes) local standard time (LST) (hours: minutes) 

Fig. 6. Comparison of aging time-scales based on number (left) and mass (right), 
showing both the time-scales due to condensation and coagulation (tn and tm) and 
the time-scales due to condensation alone ir^^^ and r^"*^). The top panels have 
the critical supersaturation separating fresh from aged particles set to 5c = 0.1%, 
while the middle panels have = 0.6%, and the bottom panels have = 1.0%. 
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Fig. 7. Transfer rates due to condensation (top) and coagulation (bottom). The 
critical supersaturation separating fresh and aged particles is set to Sc = 0.6%. The 
notation is according to Table |2j 
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Fig. 8. Day (left) and night (right) averages of the aging time-scale, as defined in 



equations (p3|)-(25). As can be see from Figure pi during the day the condensa- 



tion-induced aging dominates, so T]^°Jay is almost indistinguishable from TN^day; and 
similarly for mass. 
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